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We investigate the groundstates of weakly interacting bosons in a rotating trap as a function of 
the number of bosons, N, and the average number of vortices, Ny- We identify the filling fraction 
v = N/Ny as the parameter controlling the nature of these states. We present results indicating 
that, as a function of v, there is a zero temperature phase transition between a triangular vortex 
lattice phase, and strongly-correlated vortex liquid phases. The vortex liquid phases appear to be 
the Read-Rezayi parafermion states. 



PACS Numbers: 03.75.Fi, 73.40.Hm, 67.57.Fg 

A fundamental characteristic of condensed Bose sys- 
tems is their response to rotation p[. A transition 
to a "normal" phase might be expected at sufficiently 
high angular velocities, lu, of the container (or trap) 
by loose analogy with a superconductor in a magnetic 
field. At zero temperature this phase would constitute 
a novel uncondensed ground state. Such a regime is en- 
tered when the vortex cores start to overlap. The corre- 
sponding value of u> is unattainable with bulk 4 He, but 
may be achievable in the very dilute degenerate atomic 
gases initially explored in Ref. 0, and studied exten- 
sively in Refs. [pr|l, Apart from the identification 0] 
of the Laughlin state as the ground state at sufficiently 
high ui, work on the most interesting regime of large 
numbers of vortices has been restricted to cither mean 
field theory || or exact diagonalisation [g-g). These 
two approaches have exhibited apparently contradictory 
pictures. Within Gross-Pitaevskii (GP) mean-field the- 
ory, the groundstates are vortex lattices (distorted by 
the confinement), with broken rotational symmetry J3|. 
On the other hand, exact diagonalisations have identified 
groundstates which do not have crystalline correlations of 
vortex locations J5J; they are strongly correlated vortex- 
liquids, closely related to incompressible liquid states re- 
sponsible for the fractional quantum Hall effect |] [/]] . 

Here we present results of extensive exact diagonalisa- 
tions (EDs) that elucidate the relationship between these 
two pictures. By using a periodic geometry, we have been 
able to study systems containing many vortices up to 
boson densities far in excess of previous EDs. Our re- 
sults indicate that both crystalline and liquid phases of 
vortices exist. A clean distinction between these phases 
can only be made for a large number of vortices. In 
this limit, we argue that there is a zero-temperature 
phase transition as a function of the "filling fraction" , 
v = N/Ny, the ratio of the number of bosons, N, to the 
average number of vortices, Ny. For large ^, the ground- 
state is a vortex lattice (characterised by broken trans- 
lational/rotational symmetry). For small v the ground- 
states are strongly-correlated vortex liquids. We find that 



the vortex-liquid groundstates are related to the Read- 
Rezayi "parafermion" states |10| that were introduced in 
the context of fractional quantum Hall systems. 

In a frame of reference rotating with angular veloc- 
ity uz, the Hamiltonian for a particle of mass m in an 
(isotropic) harmonic trap of natural frequency luq is 
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The second form indicates the equivalence to the Hamil- 
tonian of a particle of charge q* experiencing an effective 
magnetic field B* = Vx(muzxr/q*) — (2mui/q*)z (the 
particle also feels a reduced xy-confinement). Of partic- 
ular importance to our discussion is the average filling 
fraction, v, for the bosons in this effective magnetic field. 
For N bosons spread over an area A, one finds 
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where Ny is the average number of vortices. For large 
Ny the vortex density is approximately uniform, and 
N v = (2muiA)/h, or (equivalently) N v = 2L/N 0, 
where L is the total angular momentum in units of Ti. 
We now introduce repulsive interactions [JL2[ 



V 



N 

9 E *(» 

»<J=1 



,). 



(2) 



with g = AttTi a/m, chosen to give the correct s-wave 
scattering length a. Throughout this work, we make use 
of the limit of weak interactions formulated in Ref. [0. 
For g <C hujQa 3 , with a the interparticle spacing, the 
bosons are restricted to single particle states in the low- 
est Landau level, and lowest oscillator state of z. For 
U) ~ Wo, the repulsive interactions give rise to the ap- 
pearance of rotating (vortex) groundstates j^-||] . 

GP theory H takes account of interactions by finding 
the fully-condensed state that minimises the total energy. 



In EDs |^|-0|, the groundstate is found by diagonalising 
(||) within the set of all states of N bosons with fixed total 
angular momentum L (L is conserved by interactions). 
An important distinction between these two approaches 
is that the GP groundstates exhibit broken rotational 
symmetry M , while the ED groundstate is necessarily an 
eigenstate of angular momentum, L. However, by per- 
forming EDs on large numbers of bosons (up to N = 30 
at L ~ 2N), we find that, as iV becomes large at fixed 
L/N , a macroscopic number of quasi-degeneracies appear 
between states with different L. This signals the emer- 
gence of broken rotational symmetry. Indeed, it appears 
from these and other |g|| studies that as N becomes 
large for fixed L/N, there is a crossover to a regime in 
which GP theory is essentially correct. We believe that 
this crossover is related to the phase transition, discussed 
in detail below, between vortex liquids at small v, and a 
vortex lattice at large v. Applying a simple Lindemann 
criterion p3| , one finds that a triangular vortex lattice 
is unstable to quantum fluctuations for v < 14. The 
crossover to GP behaviour for increasing N at fixed L/N 
is the remnant of this phase transition in a system with 
a finite number of vortices Ny — 2L/N. 

To investigate in detail the dependence of the ground- 
state on v, we have conducted extensive (Lanczos) diago- 
nalisations in a toroidal geometry [[L4J . This periodic ge- 
ometry represents the bulk of a system containing a large 
number of vortices. We consider a torus of sides a and 
b. There are then Ny = (2mujab)/h = ab(q*B* /h) vor- 
tices, which is the number of single-particle states on the 
torus in the lowest Landau level ]14J , and hence an inte- 
ger. Thus, both TV and Ny are integers, and v = N/Ny 
is a rational fraction. Finally, we classify all states by 
the Haldane momentum Jig ] , which runs over a Brillouin 
zone containing iV 2 points, where N is the greatest com- 
mon divisor of N and Ny . In the following we shall refer 
to the x and y momenta by the dimensionless vector, 
(K x ,K y ), using units of (2irh/a) and (2irTi/b). We re- 
port only positive values of K x ,K y up to the Brillouin 
zone boundary [states at (±K X , ±K y ) are degenerate by 
symmetry]. We also choose to measure energies in units 
of g/(y/t^£ 3 ), where i = v / h/(q*B*) = v / h/(2mujo) is 
the magnetic length at us — ujq. 

We start by applying Gross-Pitaevskii theory || on the 
torus. In general, the GP groundstate is a vortex lattice, 
with broken translational symmetry: the wavefunction 
is not an eigenstate of the Haldane momentum, but has 
weight at a set of reciprocal lattice vectors (RLVs). While 
the symmetry of the lattice depends, in general, on Ny 
and the aspect ratio a/b, the absolute minimum of energy 
is always obtained for a triangular vortex lattice Jig ] . 

In ED studies, the groundstate is necessarily an eigen- 
state of the Haldane momentum. The signature of trans- 
lational symmetry breaking is the development of quasi- 
degenerate levels at the set of momenta given by the 
RLVs of the broken symmetry lattice u7v . To search for 



such degeneracies, we show in Fig. [j] the evolution with 
v of the excitation energies for Ny = 8 vortices at an 
aspect ratio (a/b = y/3/4) for which the GP groundstate 
is a triangular lattice. 
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FIG. 1. Solid lines: Excitation energies at momenta mea- 
sured relative to the groundstate, for Ny = 8, a/b — \/3/4 
(inset shows the GP groundstate: dark = low boson den- 
sity). The excitation energies at the RLVs of the triangular 
lattice (filled symbols) collapse at v ~ 6, signalling the onset 
of a groundstate quasi-degeneracy; all other momenta retain 
non-zero excitation energies (two such momenta are shown as 
open symbols). Dashed Line: The excitation energy at one 
RLV (2, 0) for N v = 6 and a/b - 1/VS, showing that the 
collapse at v ~ 6 initiates an exponential decrease with v. 

A collapse of the excitation energies at the RLVs of a 
triangular lattice is observed at v ~ 6. Similar plots for 
Ny = 4, 6 indicate that the excitation energies at RLVs 
fall exponentially with v for v > 6 (shown in Fig. |lj for 
one RLV for Ny — 6). For v = 15 at Ny — 6 the exci- 
tation energies are 6 orders of magnitude smaller at the 
RLVs than at any other momentum. This strong quasi- 
degeneracy at the reciprocal lattice vectors of the lattice 
formed in GP theory indicates a strong tendency to bro- 
ken translational symmetry J17J. GP theory accurately 
describes the states at large values of v. 

We view the collapse of the excitation gaps at v ~ 6 as 
an indication, in this finite-size system, of a true phase- 
transition from translationally-invariant "vortex-liquid" 
phases, to a (triangular) vortex lattice. The phase tran- 
sition is rounded due to the finite number of vortices, 
and becomes sharper for larger Ny (over the range of Ny 
we can study). Note that we have chosen aspect ratios 
that are commensurate with a triangular lattice, which 
is likely to help stabilise the vortex lattice. Similar plots 
at other aspect ratios show transitions to a vortex lattice 
at larger values of v (up to v c ~ 15 for Ny = 4). One 
should therefore view v c ~ 6 as a lower bound on the 
critical value of v at which the transition occurs. Our 



numerical results are consistent with a transition in the 
vicinity of v c ~ 10 ± 5. 

We now turn to discuss the vortex liquids at v < 6. In 
this regime, we find incompressible liquid states similar to 
those in fractional quantum Hall systems. Some of these 
incompressible states can be accounted for by the use of a 
composite fermion construction that has previously been 
shown to describe accurately ED results on small systems 
in the disk geometry M. In the present uniform geom- 
etry, this theory predicts a sequence of incompressible 
states at v = -J^L- = ±, §, §, § . . . §, f , §, 2, oo, which 
is a bosonic version of the Jain sequence of fractional 
quantum Hall states |y| . Many of the strongest incom- 
pressible states we find cannot be accounted for in this 
way. In particular, the largest (finite) value in the com- 
posite fermion sequence is v = 2, while the transition to 
a vortex lattice does not occur until v = 6. To investi- 
gate the liquid states in this regime, we plot in Fig. || the 
energy gaps as a function of v for Ny = 6 vortices. The 
energy gap is related to the discontinuity in the chemical 
potential. To minimise finite size effects, we define the 
gap A by 
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which reduces to the standard definition as N — > oo. 
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FIG. 2. Energy gap (j3|) as a function of v for Ny = 6 vor- 
tices, at a/b = l/v3. Upward spikes signal values of v for 
which the groundstate is incompressible. The collapse of the 
gaps at v ~ 6 indicates the transition to the vortex lattice 
phase. (Inset shows the density of the GP groundstate.) 

As well as the Laughlin state at v = i B, in- 
compressible states appear clearly in Fig. H at v = 
1, |,2, |,3, |,4, |,5,6 (the loss of gaps for v > 6 is an- 
other indication of the transition to the vortex lattice, 
perhaps re-entrant around v = 6). It is not immedi- 
ately apparent how to construct incompressible states 



for this sequence of v. One possibility is that the 
vortices themselves are forming Laughlin states. This 
would provide a set of states with vortex filling fraction 
vv = \i \i \ ■ • ■ 0], and hence v = 1/Vy = 2,4, 6. . .. 
However, this construction does not account for states 
at v = 1,|,|,3, |,|,5 p0[ ]. Moreover, the trial wave- 
functions of this form that we have tested (on a disc) 
have high interaction energies: while they keep the vor- 
tices apart, they do not introduce favourable correlations 
between the bosons. States of this type are likely to de- 
scribe systems in which the underlying interactions can 
be described as repulsive two-body forces between vor- 
tices [[19 21 ll| . They do not provide an accurate descrip- 
tion in the present situation, where the interactions can- 
not be represented by pairwise vortex interactions 22 1. 

A clue to the nature of the incompressible vortex liq- 
uid states lies in the existence of an incompressible state 
at v = 1, which also cannot be accounted for in terms of 
non-interacting composite fermions. Rather, this state is 
well-described |2^| by the Moore- Read ( "Pfaffian" ) wave- 
function J24J. We find that the exact groundstate has 
large overlap with the Moore-Read state at the Haldane 
momenta for which it can be constructed on a torus p5[ . 

Motivated by this success, we have compared the in- 
compressible states at higher integer and half-integer v 
with "parafcrmion" wavefunctions introduced by Read 
and Rezayi jlQ] as generalisations of the Moore-Read 
state. These states may be represented [E6[ as a (sym- 
metrised) product of k Laughlin states via 
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where z = x + iy, and we omit the exponential factor 
of lowest Landau level states as usual. The symbol S 
indicates symmetrisation over all partitions of N parti- 
cles into sets A, B, . . . of N/k particles (we assume that 
N is divisible by k). The cases k = 1 and k — 2 corre- 
spond to the Laughlin and Moore-Read wavefunctions. 
For general k, the wavefunction (||) describes a system 
with filling fraction i/ fc ) = N 2 /2L — k/2, and is a zero 
energy eigenstate of a (k + l)-body version of the repul- 
sion (§). 

The Read-Rezayi states provide a consistent interpre- 
tation of the incompressible states in Fig. |[ they iden- 
tify the sequence of incompressible states observed in the 
EDs (y — ~ with integer k); they have large overlaps 
with the exact wavefunctions, at least up to v = 3 (the 
largest v for which we have made the comparison). We 
construct the Read-Rezayi states on the torus by diag- 
onalising the [k + l)-body force-law directly to find the 
zero energy eigenstates. In general, we find more than 
one zero energy eigenstate, and recover a total ground- 
state degeneracy on a torus of k + 1 , consistent with Ref. 



JUTf . The overlap of the exact groundstates of the two- 
body force (||) with the Read-Rezayi states are given in 
Table | for v = |. For comparison, the overlaps with the 
GP groundstate are also shown. 

In conclusion, we have shown that the groundstates 
of weakly interacting bosons in a rotating trap exhibit 
both vortex lattices and incompressible vortex liquids. 
A clear distinction between these phases appears for a 
large number of bosons, N, and vortices, Ny, and is con- 
trolled by the filling fraction v = N/Ny. Vortex liquid 
phases appear for v < v c and vortex lattices appear for 
v > v c . A Lindemann criterion suggests v c ~ 14, while 
exact diagonalisations indicate v c > 6. Current exper- 
iments |27J with N ~ 10 5 and Ny ~ 10 are deep in 
the regime in which the groundstate is a vortex lattice. 
Experiments that access the quantum-melted vortex liq- 
uid phases will require specific attention to small system 
sizes and high angular momentum. Our results indicate 
that novel correlated states emerge in this regime, which 
are well-described by the Read-Rezayi parafermion states 
whose excitations obey non-abelian statistics [Gil . 
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TABLE I. Wavefunction comparisons of the exact ground- 
states of the 2-body force-law at v — fc/2, for Nv = 6 and 
a/b = 1/V3. In each case we report: the Haldane momenta at 
which Read-Rezayi states exist, with degeneracies; the overlap 
of the exact groundstate with the Read-Rezayi state (where 
there is more than one such state, we report the total overlap 
within this set); the overlap of the exact groundstate with the 
GP groundstate (we first project the GP state onto each com- 
ponent of momentum; N/W indicates that the GP ground- 
state has no weight at this momentum). 



